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We present a theoretical study of the infrared magneto-optical properties of ferromagnetic 
(III,Mn)V semiconductors. Our analysis combines the kinetic exchange model for (III,Mn)V ferro- 
magnetism with Kubo linear response theory and Born approximation estimates for the effect of 
disorder on the valence band quasiparticles. We predict a prominent feature in the ac-Hall con- 
ductivity at a frequency that varies over the range from 200 to 400 meV, depending on Mn and 
carrier densities, and is associated with transitions between heavy-hole and light-hole bands. In its 
zero frequency limit, our Hall conductivity reduces to the fe-space Berry's phase value predicted by 
a recent theory of the anomalous Hall effect that is able to account quantitatively for experiment. 
We compute theoretical estimates for magnetic circular dichroism, Faraday rotation, and Kerr effect 
parameters as a function of Mn concentration and free carrier density. The mid-infrared response 
feature is present in each of these magneto-optical effects. 

PACS numbers: 73.20.Dx, 73.20.Mf 



I. INTRODUCTION 

Rapid progress has been achieved over the past year 
in understanding how growth and annealing conditions 
influence the properties of (III,Mn)V diluted magnetic 
semiconductor ferromagnets. These advances have led to 
the realization of samples with higher ferromagnetic tran- 
sition temperatures and conductivitiesiSi^ (III,Mn)V 
materials have normally been described using a phe- 
nomenological model^*^ in which the valence band holes 
of the host (III,V) semiconductor are coupled by ex- 
change and Coulomb interactions to Mn 2+ local-moment 
ions with spin S = 5/2. The properties predicted by this 
model are most simply understood in the strongly metal- 
lic regime for which disorder in the spatial distribution 
of the Mn 2+ ions, and other defects of the materials, can 
be treated perturbatively. These approximations lead to 
a picture of the materials in which spin-orbit coupling of 
the valence-band hole subsystem plays a key role^ in pro- 
viding detailed explanations for many qualitative effects 
discovered in experimental studies of thermodynamic and 
transport phenomena. The model can account quanti- 
tatively for the critical temperature^^ strain sensitive 
magnetic-crystalline anisotropyf£*£ anisotropic magneto- 
resistance coefficients 9 ", and the strong anomalous Hall 
effecti 10 ! 11 Golden rule estimates of quasiparticle scat- 
tering amplitudes even provide the correct order of mag- 
nitude for longitudinal dc conductivities. 9 

In this paper we discuss corresponding theoretical pre- 
dictions for the infrared magneto-optical properties of 
these materials. We evaluate magnetic circular dichro- 
ism (MCD), Faraday rotation, and Kerr effects in the in- 
frared regime for several different Mn concentrations and 
carrier densities. From the microscopic point of view, 
each of these effects reflects the non-zero value of the ac 
Hall conductivity, <J xy {ui). Our linear response theory for 
the Hall conductivity reduces in the zero frequency limit 



to the fc-space Berry's phase expression that explains dc 
Hall effect observations .ifiiAiiiSi 

In metallic ferromagnets, measurements of magneto- 
optical coefficients on band energy scales provide very 
detailed information about the influence of broken time- 
reversal symmetry on itinerant electron quasiparticle 
states. The appropriate band energy scale for the heav- 
ily p-dopcd (III,Mn)V ferromagnets, and for a number 
of other materials that have been studied recentlyi^, is 
in the infrared. For this reason, we believe that exper- 
imental infrared magneto-optical studies of (III,Mn)V 
ferromagnets are highly desirable; we expect that they 
will be carried out in the near future and that compari- 
son with the predictions presented here will be very in- 
formative in clarifying the physics of these new ferro- 
magnets. They could, for example, reveal deficiencies of 
the relatively simple theoretical formulation that we em- 
ploy. The study of the magneto-optical response of these 
ferromagnets is also potentially interesting for applica- 
tions, especially if room temperature ferromagnetism is 
achieved in the future. Magneto-optical properties of the 
closely related (II,Mn)VI diluted magnetic semiconduc- 
tor paramagnets^" have already proved useful from both 
basic science and application points of view. 

Absorption and reflection measurements in the visible 
range have been used to establish phenomenological esti- 
mates for the p-d and s-d exchange coupling constants in 
(II,Mn)VI materials, and in establishing the important 
role of valence band holes in the (III,Mn)V's. 514-1516 ' 17 
Photoemission experiments, which explore the deeper 
electronic structure have been used to explore the de- 
gree of hybridization between the underlying host va- 
lence band and Mn electronic levels, but suffer from be- 
ing surface sensitiveii2i u i2£ In the infrared regime, recent 
optical conductivity measurements have uncovered un- 
usual non-drude behavior, including an optical absorp- 
tion peak2i*22i2i24 connected to back-scattering local- 
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ization effects and to inter valence band transitions, in 
agreement with model calculations 

We organize the paper as follows. In Sec.|H]we intro- 
duce the Kubo formula description of the ac anomalous 
Hall conductivity appropriate for the studied (III,Mn)V 
ferromagnets. In Sec. IIIII we detail the model Hamil- 
tonian and approximations used in our calculations. In 
Sec. II VI we present an analytic evaluation of the anoma- 
lous Hall conductivity for the case in which disorder 
is neglected and the bands are approximated by the 
four-band spherical model. (The six-band model that 
we use for numerical calculation reduces to the four- 
band model in the limit of infinite spin-orbit coupling 
strength.) The isotropic band dispersion of this model 
makes analytic calculations possible, although they are 
still somewhat cumbersome. The details of this calcula- 
tion, which builds intuition about qualitative properties 
of the a xy (uj) curves, are relegated to an appendix. In 
Sec. [V]we present the numerical results of the full model 
Hamiltonian calculation for a xy (uj) and apply these re- 
sults to discuss all the common magneto-optical effects 
available for experiments in the present geometry. We 
summarize this work and present our conclusions in Sec. 

ED 

II. THEORETICAL APPROACH 

Our theoretical model description starts by coupling 
the host semiconductor valence band electrons, de- 
scribed within the k ■ p or Kohn-Luttinger (KL) the- 
ory, with S = 5/2 Mn local moments with a semi- 
phenomenological local exchange interaction treated at 
a mean-field level £2*2*2ii At zero temperature this gives 
rise to valence bands that are split by an effective ex- 
change field h = N Mn 2+SJ p dZ, where N Mn 2+ is the sub- 
stitutional Mn density and the strength of the exchange 
coupling is taken to be J p d = 55 meV nm -3 £L We assume 
in this paper that the magnetization is aligned along the 
growth (z) direction by an applied small external mag- 
netic field. We restrict ourselves to the T = limit, al- 
lowing us to neglect scattering off thermal fluctuations in 
the Mn moments orientation. We assume collinear mag- 
netization in the ground state, ignoring the possibility 
of disorder induced non-collinearity in the ground state 



which is known to be less likely for the strongly metallic 
(III,Mn)V ferromagnets on which we focusi^SiSS 

The linear response theory Kubo formula expression 
for the real part of the ac Hall conductivity of disorder- 
free non-interacting electrons is: 

lm[{n' k\p x \nk) {nk\p y \n' k)] 
X (uj- E r + E , r )(E ,- - E ,A ' 

V nk ' n'/e/V nk n'k> 

where \nk) are the Bloch valence-band states and E e the 
Bloch eigenenergies within k ■ p theory (we use either six 
or four band models here), m is the bare electron mass, 
/ r is the Fermi occupation number (0 or 1 at T = 0) 

for the state \nk), and p/m is the k ■ p velocity operator 
obtained 3 ^, by differentiating the k ■ p Hamiltonian with 
respect to wavevector. In the zero frequency limit Eq. ^ 
reduces to the expression used by Jungwirth et alMblL 
to explain the dc anomalous Hall conductivity of these 
materials. This recent work suggests that anomalous Hall 
effects are more quantitatively useful in characterizing 
itinerant electron ferromagnets than had previously been 
thought, at least for the present materials. In this paper 
we extend this advance to finite frequencies. 

Even though the Hall conductivity is finite in the ab- 
sence of disorder, we do anticipate that disorder will in- 
fluence the a^(oj) curves, primarily by broadening out 
features. The sources of disorder known to be relevant in 
these materials include positional randomness of the sub- 
stitutional Mn ions with charge Q = — e, random place- 
ment of interstitial Mn ions which act as double donors 
and are believed to be non-participants— in the ferro- 
magnetic order, and As anti-sites which also act as having 
charge Q = +2e and are non-magnetic. We estimate the 
influence of disorder on the valence band quasiparticles 
by calculating their lifetimes using Fermi's golden rule 
including both screened Coulomb and exchange interac- 
tions of the valence electrons with the Mn ions and the 
compensating defects^ Including disorder broadening of 
the quasiparticle spectral functions, the Kubo formula 
expression for the Hall conductivity becomes 



I 

e 2 h - 

Re[a xy (uj)} = -^r— } lm[(n' k\p x \nk) (nk\p y \n'k)] 

uum^V J— 1 

kn^n' 

x / ^f(e)A n , rk (e)Re[G^ k (e + Ku) + f(e)A n j :(e)Rc[G^: k (e - fa)], (2) 

I 

where A n g(e) = T nk -/((e — E n ^) 2 +T 2 -/4) is the disorder broadened spectral function and G rct and G adv are the 
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advanced and retarded quasiparticle Green's functions Since we are interested in the first order effects of disorder 

with finite lifetime r~i/2, obtained from the golden rule in <r xv we approximate the above expression as, 

nk y 

scattering rates from uncorrelated disorder (see Sec. II). 
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where T„ tn i = (T„ + T n >)/2 and T n are the golden rule 
scattering rates averaged over band n as in Ref. |2!| We 
use Eq. to evaluate a xy (uj) below. 
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III. MODEL HAMILTONIAN 

In the virtual crystal approximation, the interac- 
tions are replaced by their spatial averages, so that the 
Coulomb interaction vanishes and hole quasiparticlcs in- 
teract with a spatially constant kinetic-exchange field. 
The unperturbed Hamiltonian for the holes then reads 
Hq = H L + h ■ s , where Hh is the host band Hamil- 
tonian, and s is the envelope-function hole spin oper- 
ator. The host band part of the Hamiltonian is de- 
scribed via the four or six band Kohn-Luttinger model. 
Choosing the angular momentum quantization direc- 
tion to be along the z-axis, and ordering the j = 
3/2 and j = 1/2 basis functions according to the list 
(-3/2, 1/2, -1/2, -3/2; 1/2, -1/2), the Luttinger Hamil- 
tonian H L has the form 8 : 
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In the matrix (0} we have highlighted the j — 3/2 sector. 
The Kohn-Luttinger eigenenergies are measured down 
from the top of the valence band, i.e. they are hole en- 
ergies. For completeness we list the expressions which 
define the quantities that appear in H L : 



d = -^- l2 [2kl-{kl + kl)] 



(5) 



We focus here on GaAs for which 71 = 6.98, 72 = 2.06, 
73 = 2.93, and A so = 341 meV M 

We treat the effects of disorder on the hole quasipar- 
ticles through a finite lifetime scattering rate T n ^ calcu- 
lated using the Fermi's golden rule. For uncorrelated 
disorder there are two contributions to the transport 

weighted scattering rate r t = r M -" 2+ + T As -.' antl due 

. . . . . nk V> k . 

to substitutional Mn impurities and As-antisites, given 

by 
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where the scattering matrix elements are approximated 
by the expressions (in S.I. units), 

e 2 



ehosteo{\k - fc'| 2 + qj, F ) 



n 



iii, 
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^-[(7i-72)(^ + fc 2 ) + ( 7 i+ 2 72 )fc 2 
Im " 
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M n ' n i — -> -> „ ( Z nk\ Z n'k')- 

ehoBt£a[\k - fc'r + Qtf) 

Here £host is the host semiconductor dielectric constant, 
\z n ^) is the six-component envelope-function eigenspinor 
of the Hamiltonian H h , and the Thomas- Fermi screening 
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wavevector q TF = ^ g(E F )e 2 / (2e host e ) , where g(E F ) 
is the density of states at the Fermi energy, Ep. The 
inter-band scattering broadening r n! „/ in Eq. is then 
calculated by averaging r n! „/(fc) = (T n (k) + r„/(fc))/2 
over the allowed transitions between bands n and n' as 
in Ref. El 



IV. 4-BAND SPHERICAL MODEL 

In this section we briefly summarize an analytic calcu- 
lation of <j xy {u>) for a disorder free 4-band model with 
isotropic bands, the so-called spherical model. This 
model is realized by taking the spin-orbit coupling to 
infinity and taking 72 = 73 (equal to 2.5 for GaAs) in 
the Kohn-Luttinger 6-band model of Eq. 0] This yields 



H 



L-4b 



2mo 



(7i + ^72)fc 2 - 2 72 (fc-j) 5 



(6) 



with the anti-ferromagnetic coupling between the local- 
ized moments and the holes given as before by, hs z = 
(h/3)j z . From the Hamiltonian in Eq. Elone can imme- 
diately see one of the consequences of a strong spin-orbit 
coupling: for a Bloch state labeled by wavevector k, the 
spin quantization axis at h = is parallel to k. It is pos- 
sible to evaluate a x ^(cu) from Eq. [I] in this model to first 
order in h by completing a straightforward but lengthy 
exercise in degenerate perturbation theory. This calcula- 
tion is described in greater detail in appendix Here 
we simply state the final result: 



4» 



,A X y(u') 



(7) 



where the spectral function A xy (uj) is given by dif- 
ferent expressions in three different energy intervals. 
For mihE F /fi - (h/6)mih/mhh - h/2<w <mihE F /fi + 
(h/6)mi h /m h h + h/2, 
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with 



1(1 + ^) 



1(1 + f) 



(9) 



m hh = m /(7i - 272), mm = m /(7i + 272), 
£ = mih/mih, [i = m hh mih/(m h h + mih), and 
Cj = oj — hkf h /2^i, where fc;/j is the light-hole 
band Fermi wave-vector in zero exchange field. For 
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FIG. 1: Anomalous ac-Hall conductivity calculated within the 
4-band spherical model without disorder life-time broadening 
for several itinerant hole and Mn concentrations. 
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and u> = lu — hkf lh /2fj,, where k^ is the heavy-hole band 
Fermi wave- vector in zero exchange field; and A xy (u>) = 
otherwise. 

We show att(uj) for several itinerant hole and Mn con- 
centrations in Fig. ^ From the above result (and from 
the details presented in Appendix A) it is relatively sim- 
ple to see the source of the feature observed in the mid- 
infrared regime. The spectral function A xy (uS), shown 
in Fig. [21 for the parameters used in Fig. ^ nas its 
major contribution from transitions near the light-holes 
bands Fermi wave-vectors (the lower frequency peak in 
A xy {uj)) and near the heavy-holes Fermi wave- vectors 
(the higher frequency peak in A xy (u)), visible for x = 4% 
and p = 0.2nm . The transitions that contribute to first 



5 




500 1000 500 1000 

"hco [eV] tico [meV] 



FIG. 2: Spectral function A xy {uj) calculated within the 4- 
band model the itinerant hole and Mn concentrations of Fig. 
□ 



FIG. 3: Anomalous ac-Hall conductivity a xy {u) for x = 6% 
Mn concentration and p = 0.4 and 0.2nm -3 , for spherical and 
non-spherical (band-warping) models. 



order in h are between heavy and light holes with oppo- 
site polarization as shown in Appendix lAl We also note 
that there is a considerable contribution to a x ^(uj) from 
the high frequency part of the spectral function (accounts 
for rigid shifts in the low frequency range) which indicate 
the possible need to consider higher bands, maybe includ- 
ing the conduction bands, for more realistic calculations. 



V. NUMERICAL RESULTS 

The qualitative physics behind the the 4-band model 
calculation results still applies to the full model numer- 
ical calculations. However, the effects on a xy (cu) due to 
the lifetime broadening of the quasiparticles, finite spin- 
orbit coupling, and the warping of the bands (73 7^ 72) at 
higher concentrations are an important part of the quan- 
titative numerical result. As an example, Fig. [31 shows 
the anomalous ac-Hall conductivity of disordered system 
for x — 6 % and p = 0.2 and 0.4 nm~ 3 calculated using 
the 6-band model with warping (73 7^ 72) and without 
warping (73 =72), which emphasizes the importance of 
including the warping of the bands in obtaining reliable 
results which can be compared directly with experiment. 

The Hall conductivity must be non-zero in order to 
have non-zero magneto-optical effects, but most mea- 
surable quantities are also influenced by other elements 
of the conductivity tensor. The most widely studied 
magneto-optical effects are the Faraday and Kerr ef- 
fects. The Faraday effect reflects the relative difference 
between the optical absorption of right and left circu- 
larly polarized light, referred to as magnetic circular 
dichroism (MCD). In the Voigt geometry (magnetization 
aligned with axis of light propagation) and assuming a 
thin film geometry (applicable for all (III,Mn)V epilayers 



now available in the infrared regime considered here) 



MCD = 



a + - a _ lm[a xy (uj)} 
a+ + a- Re[a xx (u))]' 



(13) 



Linearly polarized light propagating through a magnetic 
medium will experience the Faraday rotation of its polar- 
ization angle and a transformation from linear to ellipti- 
cally polarized light due to MCD. The angle of rotation 
per unit length traversed, again in the thin film geometry, 
is (in cgs units) 



47T 

(1 + n)c 



(14) 



where c is the speed of light and n is the index of refrac- 
tion of the substrate, in this case GaAs with n = VlO.9. 
Perhaps the more technologically relevant magneto-optic 
phenomena is the Kerr effect, which appears in reflection 
from a magnetic medium. In this case, also within the 
Voigt geometry, the Kerr angle and ellipticity are defined 
as 



k + irjK = 



r + + T- 



(15) 



where r± are the total complex reflection amplitudes 
(with multiple scattering taken into account) for right 
and left circular polarized light. Note that the simple 
relations, 9k oc Im[t7 xy (a;)] and i]k oc Re[a xy (u>)]r^ ob- 
tained in the thick-layer limit do not apply for the typical 
thin (III,Mn)V epilayers. In Fig. 01 we show the differ- 
ent magneto-optic effects for a concentration of x — 6% 
and p = 0.4 nm -3 . The Faraday rotation in this case is 
larger than the giant Faraday rotation observed in the 
paramagnetic (II,Mn)VI's at optical frequencies^ 3 ^ and 
should be readily observable in the current highly metal- 
lic samples. The Kerr angle and ellipticity we obtain for 
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FIG. 4: Faraday and Kerr effects for x — 6% Mn concentra- 
tion and p = 0.4 nm -3 . 



(Ga,Mn)As are comparable to the Kerr effects observed 
in the optical regime in materials used for magneto- 
recording devices^ The behavior as a function of free 
carrier hole concentration can be seen in Fig. where 
the Faraday rotation angle is shown for several carrier 
concentrations. The peaks and valleys in the different 
quantities are present in all the concentrations, however 
the magnitude varies, even changing sign at several con- 
centrations and frequencies. 

Rather than presenting many different graphs for all 
the possible parameters (p,x, etc.), we direct the reader 
to a data-base located at |http : //u nixl 2 .fzu.cz /ms where 
results for these quantities, together with other physi- 
cal quantities, can be obtained and plotted vs. different 
nominal parameters^ 



VI. CONCLUSIONS 

We have presented a theory of the ac Hall effect in 
the infrared regime by extending the Berry's phase the- 
ory of the dc-anomalous Hall effect to finite frequencies 
and treating the effects of disorder through a finite life- 
time of the valence-band quasiparticles. We observe fea- 
tures (peaks and valleys) in the transverse conductivity 
in the range between 200 and 400 meV at which the con- 
ductivity changes by more than 100%. We have studied 
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FIG. 5: Faraday rotation angle for x = 6% Mn concentration 
with p = 0.2, 0.4, 0.6, and 0." 
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how these features appear in different magneto-optical ef- 
fects (MCD, Faraday rotation and Kerr effect) which are 
relatively easily measured, finding strong signals. The 
magnitude of the Faraday rotation is very large (one or- 
der of magnitude larger than that observed in param- 
agnetic (II,Mn)VI's for example) and has a nontrivial 
dependence on the free carrier concentration. The Kerr 
effect is also strong when compared to materials used 
in magneto-optic recording. The origin of the peaks is 
most easily understood within a simple 4-band spherical 
model in which transitions between heavy and light holes 
states with opposite spin-polarization give the strongest 
contribution to the anomalous transverse optical conduc- 
tivity. The four band model represents the infinite spin- 
orbit coupling strength limit of the six-band model we use 
for numerical calculations. Our use of a six-band model 
can account only for transitions within the valence band 
and not for transitions between conduction and valence 
bands. Because of this limitation, we cannot address 
the crossover between intra-band and interband contri- 
butions which are not completely separated in these ex- 
tremely heavily-doped semiconductors, something that is 
clearly desirable and should be addressed in subsequent 
theoretical work. 

Our predictions depend in intricate detail on the model 
that we have used to describe the ferromagnetism of these 
materials. The model depends most essentially on the 
assumption that the Mn impurities act as reasonably 
shallow acceptors and introduce S — 5/2 local moment 
degrees of freedom to the system. The specific calcu- 
lations presented here assume that Mn impurities and 
other scatterers in the system can be treated perturba- 
tively. This assumption enables quasiparticle scattering 
rates to be estimated in a simple way, but is a less es- 
sential part of the model. The magneto-optical proper- 
ties studied here are directly dependent on valence-band 
spin-orbit coupling, which we have argued elsewher e 4 ! 36 
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plays an essential role in understanding ferromagnetism 
in these materials. Confirmation by future experiment 
of the detailed predictions made here for the magneto- 
optical properties of these materials would further vali- 
date the approach we have taken to modeling these in- 
teresting new ferromagnets. We expect that the weak- 
quasiparticle-scattering approximations made here will 
be more reliable in more metallic samples, since the scat- 
tering rates are then smaller compared to other relevant 
energy scales, particularly the Fermi energy. We hope 
that these calculations will help motivate magneto-optic 
experiments in the infrared regime for (Ga,Mn)As and 
other (III,Mn)V ferromagnets. 
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contribution to the ac Hall conductivity calculated first 
in the exchange field within the 4-band spherical model. 
The host valence band Hamiltonian in this shown 
in Sec. IIVI is given by 
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The eigenspinors of H are given by 
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where \n) are the spinors with the axis of quantization 
along the z-direction and total angular momentum 3/2H. 
The perturbation due to the antiferromagnetic coupling 
to the localized moments is H' = hs z = (h/3)j z . The 
eigenvalues to first order in h are then given by 
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APPENDIX A: DERIVATION OF o xy {u) IN THE 
4-BAND SPHERICAL MODEL 
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We present in this appendix the details involved in where tan 29' = 2 tan 6*, hh labels heavy- holes and lh 
deriving the results shown in Eqs. I7ll2l for the anomalous labels light-holes. 



The dipole matrix elements in Eq. ^ are given by: 

(n'k\p a \nk) = ^-{z n , k \^-\z nk ) 



m(E n k - E n , k ) / d , 

~ ft K\ TirZ ) . 



y dk a 



so we can write 



lm[(n' k\p x \nk) (nk\p y \n' k)] = -p-(-E n fe - E n > k ) 2 lm 



d . . . _d_ _ 

\Zn'k\ 07 ^nkl\^n , k\ n ll %nkl 



dk x 



dk v 



where 



and similarly 



dt 



dL 



' Znk ) 



■z nk ) 



cos 6 cos 9 d 
k 09 



\Znk) 



sin 6 d 
k sin 9 36 



\znk) + cos (psm 9— \ z nk ) 



sin 6 cos 9 3 
k 39 



Znk) 



cos d) d d 
— \ Znk)+ ^6sm6-\z nk ). 



The perturbed spinor wave function can be written as 



(A5) 



(A6) 



(A7) 



(A8) 



\z nk ) =£ C "'(M)k'!> =^C: r (9,k)e- i ^ n ^ H e- i ^ h \n') = \z nk ) - -(cos^j, - sinfe)|n'), (A9) 



where j™(0) = (z n ii=kz\jz\z n k=kz) ■ Inserting Eq. IA9I into Eq. IA7I and lA"8l gives 

, 3 isind -n,^\s\ \ .cos^cosfl /. . -, ,\ , . 

z nk ) = , , . n \0z~ 3 [Q))\Znk)—i £T y[cos®j y - smdj x \\z nk ) + i\z n k) ) + cosdsi 



l 3k, 



' Znk} 



hk sin 6 
— i cos< 
hk sin t 



hk 

.„,„», x . sin </) cosy /. - . ., - >\ . , . „ O , . 

(Jz - J {0))\znk} - 1 zr [jcosdjy - smtpj x \\z„k) + i\z n k) J + smcj>sm9—\z n k), 



sm0 dk lZnk) 
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which can be inserted in Eq. IA6l to yield 



lm[{n'k\p x \nk){nk\p y \n'k)] = -j^(E nk - E n , k ) 2 {z n , k \{j z - j n (0))\z nk )lm 

+ih(z n ' k \znk) + J^\ Zn ' k ^~^^ 



cosf 



(hk) 2 sin 



{{z n 'k\ (cos cj)j y - sin 4>j x )\z nk ) 



(A10) 



where 



(z n , k \(j z - j n (0))\z nk ) = £ k)CZ(9,k){m\Q z - j"(0))|n 2 ), 



Im 



(z n > k \{cos(f>jy - sin (j>j x )\z n k) 



(^n'fcFnfc) = 2^,C ni {e,k) ^ and (z„/ fe 



88 



dk 



£<(M)- 



Here we only need to consider six transitions since we only need the n ^ n' terms and we will ignore transitions 
between bands with equal effective masses which can be shown to contribute to higher order in h. From degenerate 
perturbation theory we obtain the four eigenvectors to linear order in h: 



\k,hh±) = |k,±3/2) + ^g^|k,±l/2), 

\k,lh+) = cos0'|k,+l/2) -sin0'|k } -l/2) - h ^ sin0 [ cos fl'|k, +3/2) - sin0'|k, -3/2)], 

V3(fi./c) 2 

\k,lh-) = sin0'|k,+l/2) + cos0'|k,-l/2) - h £ am0 [ sm 0'|k, +3/2) + cos 0'|k, -3/2)], 

V3(hk) 2 

where fi = mihinhh / '{mhh — mih)- The Fermi wavevectors to first order in h/Ep for each band are given by 

h 



(All) 
(A12) 
(A13) 



k h F h± {9) = k 



hh(0) 



h 



1± cos0 andfc^ ± (6») = fc 



lh(0) 
F 



1 ± 



, , 1 cos 2 9 

6E F V 4 



(A14) 



After some lengthy algebra one obtains 



Im[(k,hh + 


\p x \k,lh+)(k,lh + 


\Py 


k,hh+)} 




i E th ~ E hh) 






lm[(k,hh + 


\p x \k,lh-)(k,lh- 


\Py 


k,hh+)} 




(E M -E+ h ) 






Im[(fc, hh — 


\p x \k,lh+)(k,lh + 


\py 


k,hh-)] 




( E hi - E hh) 






Im[(fc, hh — 


\p x \k, lh—){k, Ih - 


\py 


k,hh-)] 




( E fh - E hh) 







im 2 



( <,s 9 cos 2 9' + — - ( — sin(20) sin(20') + cos 29 cos 2 9' 
2(hk) 2 4 



cos t/ cos 



2 /)/ 



9' 3 



4 cos 29' 



3m 2 
~8/7 



cos 2 9 cos 2 9') 



hm 2 , 1 



(A15) 



■ cos v sin 



2{hk) 



2 (+- sin(20) sin(20') + cos 29 sin 2 9' 



2 n, 



cos v sin 

4 cos 29' 4 
3m 2 



cos 2 0sin^ 9'), 



(A16) 



8/i 

cos" t/ sin 



him 2 1 

2 o' + _^_( + _ sin (26») sin(20') + cos 29 sin 2 9' 
2(hk) 2 4 



2n c .\„ 2 n< 3 



cos 2 0sin^ 9'), 



4 cos 20' 4 

■- cos cos 2 0' + — — -,( — sin(20) sin(20') + cos 26* cos 2 0') 

8fi 2(hk) 2 ' 4 



(A17) 



cos tt COS 



2 /)/ 



0' 3 



4 cos 20' 



cos 2 cos 2 0') 



(A18) 
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Using Eqs. IA15IA18I we can compute directly the dc conductivity (Eq. \I\for uo = 0): 



(Txy(0) 



2e 2 h 
m 2 V 



E 

k,n>n' 



(fn'.k ~ fn,k)^[{n'k\p x \nk){nk\p y \n'k)} 



(E nk - E n , k y 



s 2 fik hh0 r 



(2nh) 4ttE f 



_ 1 pmih 



mih 



3 V JTih/> 3 m Ul + y/mihmhh 



(A19) 



in agreement with the previously derived dc-anomalous Hall conductivity^ using the Berry's phase contribution to 
the Bloch group velocity in the semi-classical equations of motion approach. 

To compute the ac-anomalous Hall conductivity given by Eq. ^ we rewrite it in terms of the spectral function 
A xy (u>) 



du 



,A xy (uj') 



with 



e 2 h x - (f n ,, k - f n ,k)Im[(n'k\p x \nk)(nk\p y \n'k)] 
A xy (uj) = -— ^ TE^. tT~~~\ 6{hhj - {E nk - E n , k )). 



(E n k ~ E n i k ) 



(A20) 



(A21) 



A xy (uj) is an odd function of lu and we need only to consider ui > 0. We need to consider three separate frequency 
ranges in what follows. First we look at the range 



mi h E F mih h h m hh E F h m hh h 

+ 7 + - < w < -, 

fi mhh 6 2 fi 6 mih 2 

and consider the different contributions to A xy (cj) from the four types of transitions, hh^ 
hh+ to lh± transitions we have 

A xy (uj\ hh-\ > I hi) 



(A22) 



Ih , separately. For 



e 2 h 



-l 



m 2 (27r 
e 2 ^2(iLj/h f 1 



1 7/ n \ ^k Im[(fc, hh + \p x \k, lh±)(k, Ih ± hh- 
,. i d(cosu) ■ 



K 2 



±- 



(2ir) 2 h 



1 cos 2 9 f cos 



d(cos 9) ( ^ cos 9 \ ^ 2 a i 



( E u - E hh) 
2 h i 



■ ? nl f COS 

8cos26>' sm 2 9' J 8 



sin 

™2 Q> 



[=F-sin(20) sin(26>') + cos 20 

ihu 4 



sin 2 9' 



sin 



2 />/ 



We can sum the two and obtain 

A xy [uj; hh-\ ► lh+) + A xy (uj; hh-\ ► Ih— ) 



e 2 v / 2^uj/h r 1 



(2ir) 2 h 



//3ft. 13 
dfcos 9) - cos 9 H [+ cos 29 + -cos 2 9 cos 2 9 
_i V 8 4ftw 8 8 

(A23) 



2/ihuj h 



(2irh) 48?r V ft 2 ftw 

For the hh— to Z/i± transition we obtain the same result, therefore within this range we have 

„2 K 



2\xhio h 
(2nh) 24tt V ft 2 feZ 



As one can see from its definition A xy {u>) changes most rapidly in the region where transitions near the Fermi 



surface are allowed. Let's next consider transitions from hh+ to lh± first in the lower range ^^-Ep 

° fi m hh 



mih h 
mhh 6 



< 



u mK E + mii^h + h 

u *" mhh 6 2 



/l />oo 
d{cos9) / dkf(9, k)5(huj — AE ± ), 
-l -% h± (0) 



(A24) 
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with 



AE* = 



± (hk) 2 hcosd h 



2[i ± 6cos26»' 2 



cos 9 and k 1 ^ (9) = kp {0} I 1 T — ^/ 1 - ^ cos 2 



The minimum of AE ± (9) + at a fixed 9 is then 

A25±(0) = ^fcr )2 ±^ C ° Sg --coBg,(A25) 
+ W 2fi F ^6cos2<9' 2 ' V ; 

where we have defined £ = — mj^///+ 1 = m/^/m/^, and 
the absolute minimum is given by 



A£5 



2/x' 



0) = :_ fc ^) 2 ± Z^_^ 

' mw, 6 2 



For hft— to lh± we have instead 



angular integration, 9, obtained by setting k = k l p so 
for hh+ to lh± 



fc - m//, h 3 „ 

no; =F W 1 cos^ i 

m hh 3V 4 



whose solution is 



• cos 9 = 



(A26) 



m^/i 6 cos 20' 2 



— cos( 



A£ d 



ff_,ih(v) 2 ± ma _ h 
2fi F m hh 6 2 



cos#± 



to±|^//i 2 (l + f ) - 3fc 2 d; 2 
1(1 + 1) 



1 - A^A27) 



Let frw = ^k l p {0) + Co where frw will be of the order of A similar procedure for the transitions from hh- to lh± 
h. For an Co that is too small there will be a limit on the yield cos 9± = — 1 + A T . 

I 

Combining the contributions for each transition we then obtain for mihEp/^i — (h/6)mih/mhh — h/2 < lo < 
mi h E F /fi+ (h/6)mi h /m hh + h/2 



A xy (to) = 



/ d(cos 9)A xy (u>, cos 9; hh^ >lh+)+ / d(cos 9)A xy (lo, cos 9; hh^ ► Ih— ) 

J1-A+ Jl-A- 



-i+A_ 



d(cos 9)A xy (oj, cos 9; hh ► + 



-1+A+ 



d(cos 9) A xy (co, cos 9; hh ► i/i— ) 



l v / 2 f uo/h f 1 



{2ir) 2 h 



f d( cos 9) (- cos 9 cos 2 9 

Jl-A + V4 



2a ' 1 — [ — sin(26>) sin(2(9') + 2 cos 2(9 cos 2 6" 
4nu> 2 



cos t/ cos 



4 cos 26*' 



2 /)/ 



■ cos 2 (9 cos 2 9'} 



e 2 v /2/iw/fc Z 11 



(27r) 2 fc 



/" /3 hi 

/ d(cos6») -cos6>sin 2 (9' + [+- sin(26») sin(2(9') + 2 cos 2(9 sin 2 6" 

Ji-A- V 4 4ficj 2 



cos 2 9 sin 2 6" 3 2 n . , _ 

— cos 2 9 sm 2 9' 

4cos26>' 4 



e 2 yj2^,uo/h 
(2nj 2 h~ 



3 ? h r 7 r, „ , 
8" + 4^ [ (T - 2M] 



l-A , 



1 - -u 2 



— -| arcsin(v / 3u/2) 



l-A, 



1-A_ 



Ahuj 



7u 3 
12" 



1 3 2 
1 li 

4 







1-A+- 




13u 2 \\ 












1-A__ 



A similar procedure for the upper range h/6+ (h/2)rrihh Imih yields A xy (u>) given in Eg. 1111 For 

m hh E F /mu - h/6 - (h/2)m hh /mi h <Lo<m hh E F /mu + 
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any other value of to A xy (u>) = 0. 
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